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HIGHLIGHTS 


• A comprehensive model is proposed for the Lurgi fixed-bed gasifier. 

• A strategy is proposed to accelerate the convergence of the model solutions. 

• Satisfactory agreement between model predictions and industrial data is obtained. 

• An exergy analysis is applied to both the Lurgi gasifier and the gasification system. 
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This paper presents a comprehensive steady state kinetic model of a commercial-scale pressurized Lurgi 
fixed-bed dry bottom coal gasifier. The model is developed using the simulator Aspen Plus. Five sequen¬ 
tial modules: drying zone, pyrolysis zone, gasification zone, combustion zone and overall heat recovery 
unit, are considered in the main process model. A non-linear programming (NLP) model is employed to 
estimate the pyrolysis products, which include char, coal gas and high-weight hydrocarbons/distillable 
liquids (tar, phenol, naphtha and oil). To accelerate solution convergence, an external FORTRAN subrou¬ 
tine is used to simulate the kinetics of the combustion and gasification processes which are formulated in 
terms of a series of continuous stirred-tank reactors. The model is validated with industrial data. The 
effects of two key operating parameters, namely oxygen/coal mass ratio and steam/coal mass ratio, on 
the thermodynamic efficiencies of the Lurgi gasifier and the gasification system as a whole are investi¬ 
gated via extensive simulation studies. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 

Fixed-bed gasifiers, which offer a high degree of feedstock flex¬ 
ibility (coal, petroleum coke, biomass, municipal waste, etc.), are 
the most widely used technology for producing chemicals, power, 
clean alternative fuels, etc. Among commercial fixed-bed gasifiers, 
the typical and recognizable technology is crushed coal pressurized 
gasification, namely Lurgi gasification, which can handle low rank 
coals with high-ash, high-moisture and stable chemical bonds (e.g., 
lignite, sub-bituminous and bituminous coals) [1]. The Lurgi gas¬ 
ifier is relatively simple, reliable, highly efficient, and can produce 
a gaseous fuel with relatively high CH 4 concentration and high ra¬ 
tio of H 2 /CO. The Lurgi gasification technology has attracted signif¬ 
icant commercial interest in China owing to its suitability for 
production of synthetic natural gas (SNG), pipeline gas, methanol 
and F-T synthetic oil. 


* Corresponding author. Tel.: +86 18611446202. 
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A comprehensive mathematical model of fixed-bed gasifier can 
provide an effective engineering tool for thermodynamic assess¬ 
ment, abnormal state diagnosis and operation optimization. Since 
the 1970s, several attempts have been made to develop fixed- 
bed gasifier models that are largely based on mechanistic consider¬ 
ations. Examples include several two-dimensional models pro¬ 
posed by Yu et al. [2] Yoon et al. [3,4] and Hobbs et al. [5,6] a 
FORTRAN-based kinetic model developed by Wen et al. [7] a heter¬ 
ogeneous model proposed by Joseph and co-workers [8,9] and a 
two level (particle level and reactor level) model developed by Gri- 
eco and Baldi [10]. In these mechanistic models a large number of 
mathematical equations are used to describe complex momentum 
flows, mass flows, physical properties and chemical kinetics. These 
models have been used to predict the axial profiles of gas and solid 
phase temperatures, flow rates and product distributions. The 
more recent mechanistic modeling approach is geared toward 
multi-dimensional models that incorporate an overall kinetic 
model of pyrolysis/gasification/combustion processes under differ¬ 
ent operational and control conditions. However, these rigorous 
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Nomenclature 


a atom 

A area (m 2 ) 

At all atom set in pyrolysis reaction 

Q molar concentration of component i (mol m -3 ) 

d diameter (cm) 

E energy/exergy (MJ kg -1 ) 

G slack variable among coal gas, coal liquids and coal 

water 

H enthalpy (kj kg -1 ) 

h gasifier height (m) 

HHV high heating value (MJ kg -1 ) 

/ irreversibility (kj kg -1 ) 

/< coefficient or constant 

L slack variable of coal liquids 

LHV low heating value (MJ kg -1 ) 

mf mass flow (kg h _1 ) 

m RCSTRs’ numbers in combustion process 

n(C) mole fraction of carbon in fuel (%) 

n(H) mole fraction of hydrogen in fuel (%) 

n(O) mole fraction of oxygen in fuel (%) 

ME mean error 

MW molar weight (g mol -1 ) 

N molar flow rate (mol h _1 ) 

n RCSTRs’ number in combustion process 

Na number of atoms 

O/C oxygen to coal mass ratio 

P pressure (atm) 

P t total pressure (atm) 

Q heat (kj (kg h) _1 ) 

R slack variable of coal gas 

RE relative error 

RCSTR continuous stirred-tank reactor 
s species 

S entropy (kj (kg K) 1 ) 

S/C steam to coal mass ratio 

T temperature (K) 

T dp average temperature of drying-pyrolysis zone (K) 

U heat transfer coefficient (J (m 2 h K) -1 ) 

V volume (1<N m 3 ) 

Vp vapor-phase set in pyrolysis process 
w mass fraction (%) 

w* mass portion of volatile product at studied pressure 
w* r mass portion of volatile product at very low pressure 
(P< 0.001 atm) 


w** mass portion of volatile product at very high pressure 
(P> 100 atm) 

X conversion ratio (%) 

y volume fraction (%) 

y ratio of radius of unreacted core (cm) 

Greek letters 


a 

pressure constant 

P 

correlation between LHV and exergy 

£ 

porosity of ash 

V 

weight factor of L 

y 

weight factor of G 

K 

weight factor of R 

2 

a CO/C0 2 -related term 

rj 

efficiency (%) 

Subscripts 

0 

ambient conditions (298.15 K, 1 atm) 

CG 

cold gas efficiency (%) 

dash 

ash diffusion constant 

diff 

gas film diffusion 

dp 

drying-pyrolysis zones 

g 

gas phase 

jacket 

water jacket of gasifier 

particle 

char particle 

pyr 

pyrolysis reaction 

Pyr,i 

atom i in the pyrolysis reaction 

Pyrj 

gaseous species j in the pyrolysis reaction 

s 

chemical reaction constant 

tr 

heat transfer (kj (kg h)" 1 ) 

uc 

unreacted coal 

wall 

exterior wall of gasifier 

wgs 

water gas shift 

X 

exergy (kj kg -1 ) 

Superscripts 

ch 

chemical exergy 

com 

combustion zone 

dp 

drying-pyrolysis zones 

gasi 

gasification zone 

ph 

physical exergy 


models containing partial differential equations often require a 
numerical solution. In addition, the utility of these models is some¬ 
what limited because they require inordinate effort to generate 
substantial data for model validation. 

One possibility to reduce the complexity of the mechanistic 
modeling approach is to develop the simplest possible model that 
incorporates the principal gasification reactions and gross physical 
characteristics of the reactor. Such an approach may be conve¬ 
niently realized using the process simulator Aspen Plus. Aspen Plus 
has been widely employed in the simulation of carbonaceous fuel 
conversion processes. For example, it has been successfully used 
to simulate an IGCC (integrated gasification combined cycle) plant 
[11-13] a poly-generation system based on Texaco/Shell/Siemens 
gasification technology [14-16], a coal-to-methanol process [17], 
a direct coal liquefaction process [18], a compartmented fluid- 
ized-bed coal gasifier [19], a pressurized ash agglomerating fluid- 
ized-bed gasifier [20], and a coal hydro-gasification process [21]. 
It is considered to be an excellent process design and simulation 
tool because of its ability to provide comprehensive property data 


and simulate a variety of steady-state processes involving many 
interconnected units. However, the vast majority of the aforemen¬ 
tioned Aspen Plus gasifier models focus on the entrained gasifier or 
fluidized-bed gasifier type which can be represented by the specific 
chemical equilibrium-based Gibbs reactor. In a typical fixed-bed 
gasifier, owing to the relatively low temperature in the upper sec¬ 
tion of the gasifier where the coal is dried and pyrolyzed (devola¬ 
tilized), the pyrolysis products, including coal gas and high-weight 
hydrocarbons (tar, phenol, oil, etc.), are difficult to crack or burn. As 
a result, the pyrolysis products mix with the upflow gas and flow to 
the downstream units where they are recovered as by-products. 
Therefore, the yields of the pyrolysis products will affect the com¬ 
position of the syngas. This pyrolysis phenomenon is often ignored 
by most of the commonly used Aspen Plus fixed-bed gasifier mod¬ 
els. To date, very few attempts to construct a detailed pyrolysis 
model have been made due to a lack of pertinent data [5]. 

Further, there are no standard functions available in Aspen Plus 
for simulating the gas-solid countercurrent contacting pattern 
within a fixed-bed reactor. Table 1 gives an overview of some As- 
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Table 1 

Overview of Aspen Plus fixed-bed gasifier models. 3 


Model name 

Feedstock 

Developers/year 

Ash 

Unique features 

Ref. 

RGAS 

Lignite coal 

Rinard and Benjamin/1985 

Dry ash 

Countercurrent flow; kinetic routines available 

[22] 

WEN 

Bituminous/sub-bituminous/lignite coal 

Wen et al./1984 

Dry ash 

Kinetic model 

[7] 

MIT 

Bituminous coal 

MIT Technology/1985 

Slagging 

Nonkinetic model; limestone fluxing agent 

[23] 

ORNL 

Bituminous coal 

Oak Ridge National Laboratory/1985 

Dry ash 

Nonkinetic model 

[23] 

Aspentech 

Lignite coal 

Aspentech/2010 

Dry ash 

Countercurrent flow; kinetic routines available 

[24] 

BDCL 

Lignite coal/dried Poplar 

Li et al./2010 

Dry ash 

Countercurrent flow; nonkinetic model 

[25,26] 

MSW 

Municipal solid waste 

Chen et al./2013 

Residue 

Nonkinetic model 

[27] 


3 Adapted from Ref. [23]. 


pen Plus fixed-bed gasifier models. Li et al. [25,26] developed an 
Aspen Plus simulation model of a syngas chemical looping (SCL) 
gasification process using a multistage equilibrium model which 
showed good agreement with the characteristics of chemical reac¬ 
tions and mass transfer in the fixed-bed gasifier. However, the 
intricacies of the physical and chemical phenomena occurring in¬ 
side the fixed-bed gasifier as well as the gas-solid countercurrent 
movement make the kinetic simulation of such process a difficult 
task. To simulate the kinetic behavior of the Lurgi fixed-bed gas¬ 
ifier, a multistage kinetic model formulated in term of a series of 
continuous stirred-tank reactors (RCSTRs) may be used [24]. 
Although this approach ignores the temperature difference be¬ 
tween the gas and solid phases in the actual gasifier, it is still pos¬ 
sible to obtain reasonable results that reflect the basic chemical 
reactions and mass transfer process. However, existing fixed-bed 
kinetic models based on Aspen Plus suffer from the drawbacks of 
time-consuming calculation and non-convergence of the model 
solutions [23]. 

Recently, AspenTech, the developer of Aspen Plus, has devel¬ 
oped a kinetic model for moving bed coal gasifier [24]. However, 
its practical utility is somewhat limited. The most obvious limita¬ 
tion is that the moving bed coal gasifier model exhibits conver¬ 
gence difficulties when applied to different numbers of RCSTRs. 
The simulation results of the model [24] deviate from data taken 
from the literature [7], and the simulated CO concentration is 
noticeably higher than that of a practical fixed-bed coal gasifier 
(10-20%). In addition, the moving bed coal gasifier model uses 
experimental data from the work of Suuberg et al. [43] to deter¬ 
mine the product distribution during pyrolysis. This approximate 
approach fails to consider the overall atom balance, pressure effect 
and some vapor-phase products, e.g., HCN, COS, NH 3 , and the high- 
weight hydrocarbons (tar, phenol, naphtha and oil). Note that the 
high-weight hydrocarbons are considered important by-products 
in a coal-to-SNG plant. This paper presents a novel method to ad¬ 
dress these deficiencies simultaneously, and this is achieved by 
integrating GAMS modeling results and external FORTRAN subrou¬ 
tine into the Aspen Plus simulation framework. 

A related issue is the thermochemical conversion efficiency of a 
gasification process. Prins et al. [28,29], Rao et al. [30] and Jurascik 
et al. [31] have evaluated the performance of gasification systems 
in terms of different thermodynamic efficiency indicators. An 
effective tool to estimate the efficiency of a gasification process 
and determine the energy quality and usefulness is exergy analy¬ 
sis. This thermodynamic analysis allows one to determine the 
maximum performance of a system and the sources of thermody¬ 
namic irreversibility. Previous thermodynamic analyses have 
mostly focused on the effect of operating parameters such as 
equivalence ratio, steam injection rate, gasification temperature 
and pressure on the solid fuels gasifier [20,32]. However, the gas¬ 
ifier is just one part of the gasification system which consists of 
other important components such as steam, oxygen and electricity 
generation units. It has been shown that the exergy efficiency of a 
fluidized bed gasifier was different from that of the gasification 


system as a whole [28,29]. As such, the exergy analysis reported 
here evaluates the effect of different operating conditions on the 
efficiencies of a fixed bed gasifier as well as the gasification system 
as a whole. 

This work presents a practical mathematical framework to 
model, simulate and evaluate the performance of a commercial- 
scale Lurgi fixed-bed coal gasifier. The major novelties of this con¬ 
tribution are summarized as follows: 

• Based on experimental data, we formulate a non-linear pro¬ 
gramming (NLP) mathematical model that accounts for pres¬ 
sure effect to estimate the actual yields of pyrolysis products, 
including char, coal gas and high-weight hydrocarbons/distilla¬ 
ble liquids (tar, phenol, naphtha and oil). 

• To accelerate model convergence, this work presents a solution 
strategy to reduce the computational time and scale and com¬ 
plexity of the algorithm. In the simulation process, different 
external FORTRAN subroutine codes are used to link m number 
of RCSTRs in the char combustion process and n number of 
RCSTRs in the gasification process. 

• The exergy loss of each process component in the Lurgi gasifier, 
such as pyrolysis, combustion and gasification as well as the 
exergetic efficiency of the overall gasification system are 
determined. 

The remainder of this paper is organized as follows. The Lurgi 
fixed-bed gasifier under study is described in some detail in Sec¬ 
tion 2. Section 3 presents the model formulation, including the coal 
pyrolysis process model, models of char gasification and combus¬ 
tion processes, and heat transfer model. Section 4 introduces four 
kinds of exergetic efficiency indicators. In Section 5, we report 
the modeling results obtained along with some discussion of ther¬ 
modynamic analysis. Finally, the conclusions of this work are sum¬ 
marized in Section 6. 


2. Description of the Lurgi fixed-bed gasifier 

Due to strong market interest, the energy-chemical industries 
in China are investing heavily in clean coal conversion processes, 
especially the development of commercial coal-to-SNG demon¬ 
stration projects [33]. The crushed coal pressurized Lurgi gasifier 
has been successfully tested in a 4 billion N m 3 SNG/year coal-to- 
SNG plant located in Liaoning province. Fig. 1 depicts a schematic 
drawing of this updraft pressurized Lurgi gasifier analyzed in this 
paper. The Lurgi gasifier in the SNG plant consists of nine parts 
(coal lock hopper, coal chute, water jacket, gas-liquid separator, 
scrubbing cooler, coal distributor, coal gasifier bed, ash grate, and 
ash lock hopper) and has a diameter of 3.85 m, a bed height of 
12.5 m and a production capacity of 40 km 3 SNG/h (dry basis). As 
illustrated in Fig. 1, the whole gasification process in the gasifier 
bed can be divided into five individual zones: drying, pyrolysis, 
gasification, combustion and ash. 
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Fig. 1. Schematic diagram of pressurized Lurgi gasifier. 


In this gasifier, the solid phase (coal) is fed at the top of the gas¬ 
ifier lock hopper and moves downward under gravity flow coun¬ 
tercurrent to the upward gaseous stream. The first zone is the 
drying process, in which the total physical moisture bound in the 
coal is released into the gaseous phase, and the coal is pre-heated 
to 473-773 K. Further heating at higher temperature takes place in 
the pyrolysis zone where the dry coal undergoes pyrolysis/devola¬ 
tilization in the absence of oxygen to produce volatiles consisting 
of a mixture of char, coal gas and liquids, as depicted in Reaction 
Rl. 

Dry coal pyr ^!> sls char + coal gas + liquids (Rl) (1) 

Char is defined in this study as un-distillable material which 
only consists of fixed carbon. Coal gas refers to those components 
lighter than C 6 , including CO, H 2 , H 2 0, CH 4 , C0 2 and C n H m (n ^ 6). 
The term liquids represents distillable organic liquids (high-weight 
hydrocarbons), including tar, phenol, naphtha and oil, whose 


molecular weights are larger than that of C 6 . Note that the liquids 
components have indefinable compositions and liquid-like proper¬ 
ties at ambient conditions. 

After the drying and pyrolysis zones, the formed char flows to 
the gasification zone where it reacts with high temperature gas¬ 
ifying agents (steam, hydrogen and carbon dioxide), see reactions 
R2-R4 in Table 2. In addition, the CH 4 shift (R6) and CO shift (R7) 
reactions are also taken into account in the present gasification 
model. Char gasification involves a series of endothermic reac¬ 
tions which yield combustible gases such as CO, H 2 and CH 4 . 
After gasification, the remaining char flows to the combustion 
zone where combustion reactions (R5 and R8) and CO shift reac¬ 
tion (R7) occur. Char partial oxidation (R5) is a dominant reaction 
taking place in the combustion zone, which generates almost all 
the heat needed by the endothermic reactions. Finally the uncon¬ 
vertible ash and unreacted char fall into the ash zone. It is as¬ 
sumed that negligible reactions take place in the ash zone. The 
reactions R2-R8 listed in Table 2 are considered in the gasifica¬ 
tion and combustion zones. Their chemical reaction rate expres¬ 
sions and corresponding reaction rate constants are given in 
Section 3.2. 


3. Model formulation 

The process simulator Aspen Plus™ (Version 7.3) is used to sim¬ 
ulate the Lurgi gasifier and the whole gasification system. Addi¬ 
tionally, FORTRAN 11.0 and GAMS 22.3 are used to facilitate the 
modeling of the chemical kinetic reactions and coal pyrolysis pro¬ 
cess, respectively. Table 3 lists the properties of the coal examined 
in this study. It is a highly volatile lignite produced by ShengLi 
Coalfield, Inner Mongolia, China. 

In Aspen Plus coal can be treated as a mixture composed of a 
series of stable elementary substances, such as C (s) , H, N, Cl, S (s ), 
O and ash forming elements because there is no molecular formula 
for describing the coal chemical composition due to its natural 
structural complexity and composition diversity [35]. In our mod¬ 
el, the MIXCINC stream class is chosen as the global stream class. It 
is well recognized that coal and ash are treated as non-conven- 
tional components in the Aspen Physical Property Database. The 
two built-in models, HCOALGEN and DCOALIGT, are used to calcu¬ 
late the enthalpy and density of non-conventional components, 
respectively. As for the MIXED conventional components and 
CISOLID components, their physical properties are calculated by 
the PR-BM physical method dynamically as functions of tempera¬ 
ture and pressure. Fig. 2 shows the Aspen Plus simulation model 
and the calculation procedure used. According to the actual pro¬ 
cesses, five Aspen Plus blocks are used to simulate the gasifier-dry¬ 
ing zone, pyrolysis zone, gasification zone, combustion zone and 
overall heat recovery unit, in which the heat recovery is achieved 
by means of a water jacket and a heat recovery steam generator 


Table 2 

Reactions considered in the gasification and combustion zones. 


Reaction 

Reaction heat 3 , MJ mol 1 

Reaction name 

Reaction no. 

Heterogeneous reaction 

C + H 2 0 -> CO + h 2 

+131 

Char gasification 

R2 

C + 2H 2 -> CH 4 

-75 

Methanation 

R3 

C + C0 2 -> 2CO 

+172 

Boudouard 

R4 

C + x° 2 ^(x- 1)C0 2 + (1-1)CO, X e [l,2] b 

- 

Char partial oxidation 

R5 

Homogeneous reaction 

CH 4 + H 2 0 -> CO + 3H 2 

+206 

CH 4 shift 

R6 

CO + H 2 0 -► C0 2 + h 2 

-41 

CO shift 

R7 

H 2 + 0.5O 2 -> H 2 0 

-242 

H 2 combustion 

R8 


a Negative sign denotes an exothermic reaction and positive sign denotes an endothermic reaction. 

b C0/C0 2 = Ae~ E / RT = 2(2 — l)/(2 — 2), where A = 10 5 , E = 38 kj mol -1 at low pressure, and A = 10 35 , E = 50 kj mol -1 at high pressure [34]. 
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Table 3 

Properties of ShengLi coal. 


Proximate analysis, wt% 



Moisture, wet 

17.81 

Fixed carbon, 
dry 

49.5 

Volatile matter, 

35.88 

Ash, dry 

14.62 

dry 




Ultimate analysis, wt% 



Ash, dry 

14.62 

C, dry 

63.55 

H, dry 

3.89 

N, dry 

0.74 

Cl, dry 

0 

S, dry 

1.9 

0, dry 

15.3 



Particle size distribution, % 



0-10 mm 

5 

10-25 mm 

70 

25-50 mm 

20 

50-100 mm 

5 

Molecularity, daf/ 

C 6 .i 97 H 4 . 524 No.O 6 l 9 S 0 .O 69 O 1 .123 

HHV ad , Mj kg" 1 

19.87 

d 




Ash fusibility a , I< 

DT 

1500 

ST 

1510 

HT 

1510 

FT 

1530 


a DT: deformation temperature; ST: softening temperature; HT: hemispherical 
temperature; FT: flowing temperature. 


(HRSG). A description of the different block conditions and opera¬ 
tions in the Aspen Plus model is given in Table 4. 

To simplify the simulation of the Lurgi gasifier, the following 
basic assumptions are adopted in this paper. 

• The combustor and the gasifier operate under a steady state, 
and mixing of gas-solid in each reaction zone is perfect. 

• Ash in coal as well as in bed particles (sand) is inert and does 
not participate in any chemical reactions. 


• The net thermal effect of the pyrolysis process is neutral. No 
high-weight hydrocarbons are recycled to the gasifier, and the 
char contains only fixed carbon (graphite). 

• The gas products of coal gasification are H 2 , CO, C0 2 , CH 4 , H 2 0, 
C 2 H 6 , N 2 , H 2 S, NH 3 , HCN, COS, Ci 4 H 10 , C 6 H 6 0, C 7 H 8 and C 7 H 10 . 
The last 10 products can only be generated from the pyrolysis 
process. 

• The pressure loss in the gasifier is not taken into account in the 
simulation. 

• All char particles are assumed to have the same spherical 
dimensions. There is no temperature gradient between the gas 
and solid phases in the radial direction of the bed. 

3.2. Coal pyrolysis process model 

It is well known that the “RYIELD/DECOMPOSITION” block in 
Aspen Plus is usually used to model the pyrolysis/devolatilization 
process in the Shell and Texaco gasifiers, in which the coal is con¬ 
verted into a mixture composed of a series of stable simple sub¬ 
stances, including C (s ), H 2 , N 2 , Cl 2 , S( s) , 0 2 , H 2 0 and ash forming 
elements [20,36-39]. Additionally, some studies have developed 
an equivalent pyrolysis model that uses correlation equations for 
calculating the gaseous volatile matter (CH 4 , H 2 , C0 2 , CO, H 2 0, 
and tar) in a fluidized bed gasifier [40,41]. 

Several correlations derived from typical data for different coal 
types have been developed to predict the vapor-phase products of 
the pyrolysis step in a fixed bed gasifier model [42]. These correl¬ 
ative methods do not require detailed data about ultimate analysis 
and pyrolysis products of different coal types but they will not con¬ 
sistently close the overall atomic balance and mass balance during 
calculation. As such, the predicted results of the aforementioned 




Fig. 2. (a) Aspen Plus simulation calculation procedure and (b) Aspen Plus simulation model of pressurized Lurgi gasifier. 
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Table 4 

Description of block conditions and operations in the Aspen Plus model. 


Block ID 

Block type 

Description 

DRYING 

RYield 

(1) Reduces moisture content of coal 

(2) Temperature = 473 K 

(3) Pressure = 40 atm/35 atm a 

PYROLYS 

RYield with built-in FORTRAN code 

(1) RYield reactor converts non-conventional coal into conventional components 

(2) Temperature = 773 K 

(3) Pressure = 40 atm/35 atm a 

GASIF 

HERARCHY with external FORTRAN code 

(1) Uses RCSTRs in series to represent the char gasification reaction 

(2) Pressure = 40 atm/35 atm a 

COMB 

HERARCHY with external FORTRAN code 

(1) Uses RCSTRs in series to represent the char combustion reaction 

(2) Pressure = 40 atm/35 atm a 

SEP 

SEPRATOR 

(1) Separates water/gas from coal/char by specifying split fractions 

(2) Pressure = 40 atm/35 atm a 

(3) SEP1/2 temperature = 473 K/773 K 

MIX 

Mixer 

(1) Pressure = 40 atm/35 atm a 
(1) Heat losses = 0 

HRSG 

HeatX 

(1) Recovers waste heat from crude syngas 

(2) Outlet syngas temperature = 473 K 

H-JACKET 

Heater 

(1) Recovers heat from furnace wall and generates medium-pressure saturated steam for gasification 

(2) Temperature = 523 K 

(3) Pressure = 35 atm 


a 40 atm is for Runs 1 and 2, 35 atm is for Runs 3 and 4 (the operation conditions of Runs 1-4 are given in Section 5.1). 


correlative approach will show deviation from observed results. 
Compared with other types of gasifiers, the gas-phase species 
and solid-phase char in the Lurgi gasifier are in countercurrent 
movement as well as in a relatively low operational temperature 
range, so the high-weight hydrocarbons to be decomposed or com¬ 
busted before being entrained out of the gasifier by the upward 
syngas are rather different. Generally, the high-weight hydrocar¬ 
bons include tar, phenol, naphtha and oil, accounting for about 
5% of the total syngas product weight. Thus, the pyrolysis step is 
one of the most important steps in fixed-bed gasification, which 
will exert a significant influence on the final syngas product 
composition. 

The coal pyrolysis phenomenon is still not completely under¬ 
stood. It is difficult to relate the reaction mechanism to fuel prop¬ 
erties and such operating conditions as temperature, pressure, 
particle size, constituents of carrier gas, residence time as well as 
heating rate. The fundamental thermodynamic/dynamic data are 
rather scarce and difficult to obtain. Loison and Chanvin [42] pro¬ 
posed some correlations to estimate mass fractions of the most 
important gaseous products in coal pyrolysis. The correlations re¬ 
quire only an approximate analysis of coal without accounting 
for the effects of some of the operating parameters noted above. 
Suuberg et al. [43] correlated experimental results of lignite pyro¬ 
lysis at ambient pressure using a complex model that considers the 
distribution of activation energies. Pickett [44] provided some cor¬ 
relations that characterize the mass relationships among the tar, 
phenol, naphtha and oil components but the empirical equations 
only consider the mass of feed coal. Until now, published studies 
on the modeling of the pyrolysis process using Aspen Plus or other 
commercial software have not attempted to model the volatile 
product composition. 

Yoon et al. [3,4], Adanez and Labiano [45] and Baliban et al. [46] 
concluded that using models that are based on pyrolysis experi¬ 
ments is simpler and more practical in predicting the actual prod¬ 
ucts of coal pyrolysis. Coal devolatilization data could reflect the 
distributions of gaseous products and distillable liquids in coal 
pyrolysis. Incorporating the knowledge of coal devolatilization in 
the pyrolysis process model is crucial for obtaining a precise pre¬ 
diction [47]. In this paper, the amount of each pyrolysis product 
is based on literature data [24,42,43], as shown in Table 5. The ulti¬ 
mate analysis and proportion of high-weight hydrocarbons/distil- 


Table 5 

Typical coal pyrolysis/devolatilization data. 


Fractional distribution of volatiles 


Ref. 

Distribution of coal gas, wt% 

55.47 

[42,43] 

CO 

20.95 


co 2 

24.77 


ch 4 

3.82 


h 2 

1.48 


C 2 H 6 a 

4.45 


n 2 

- 


h 2 s 

- 


cos 

- 


Tar, oil, phenol, naphtha, wt% 

15.91 

[4,42] 

Water, wt% 

28.61 

[5,42,43] 


a C 2 H 6 denotes the total fraction of hydrocarbons whose molar weights are 
higher than C 2) e.g., C 2 H 4 , C 2 H 6 , C 3 H 8) C 4 H 10 , etc. 


lable liquids were obtained from actual laboratory-scale plants. 
As listed in Table 6, the tar, phenol, naphtha and oil are respec¬ 
tively represented by Ci 4 Hi 0 (pyrene), C 6 H 6 0 (phenol), C 7 H 8 (tolu¬ 
ene) and C 7 H 10 (2-norbornene), since these conventional 
components are available in the Aspen Plus property database 
and have physical properties such as C/H or molar weights that 
are similar to the corresponding high-weight hydrocarbons. 

The literature data listed in Table 5 have been obtained from 
different experimental studies. Their operating conditions (e.g., 
gasifier pressure and average temperature) are different from those 
of the fixed-bed gasifier under study here. As mentioned above, it 
is difficult to close the overall atomic balance and mass balance 
simultaneously during calculation if we use these data directly. 
However, using an approximate mathematical model for parame¬ 
ter estimation based on these experimental data will give an in¬ 
sight into the initial composition of high-weight hydrocarbons 
and provide reasonably accurate effluent flow rates for the final va¬ 
por-phase products. 

Coal pyrolysis is assumed to be instantaneous without thermal 
effects and to occur during heating. In a fixed-bed gasifier, the coal 
undergoes a slow heating pyrolysis process (heating rate <5 K s -1 ) 
and the bulk reaction temperature (700-900 K) does not change 
significantly under normal operating conditions. Most experimen- 
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Table 6 

Ultimate analysis and proportion of distillable liquids. 



Tar 

Phenol 

Naphtha 

Oil 

Ash, wt% 

1.07 

Trace 

Trace 

0.026 

C, wt% 

91.55 

76.58 

87.817 

87.511 

H, wt% 

4.95 

6.42 

11.04 

10.576 

N, wt% 

1.63 

0 

0.098 

0.899 

Cl, wt% 

0 

0 

0 

0 

S, wt% 

0.80 

0 

1.044 

1.013 

O, wt% 

0 

17.00 

0 

0 

Molecularity, daf 

C 7 . 62 H 4 . 90 N 0 . 11 S 0.03 

C6.38H6.360i.06 

C 7 . 56 H 7 . 93 N 0 . 01 S 0.03 

C7.29H 10.47N 0 . 06 S 0.03 

Density, g cm -3 

1.02 

1.063 

0.803 

0.903 

Freezing point, K 

310.1 

312.03 

242.6 

245.1 

Fraction, wt% 

33.67 

12.79 

11.94 

41.60 

Alternative ID 

C 14 H 10 

c 6 h 6 o 

c 7 h 8 

C 7 H 10 


tal data on volatile yields were measured at 1 atm in the laboratory 
while the pressure condition in practical fixed-bed gasifiers is typ¬ 
ically higher than 20 atm. The pressurized condition will influence 
the diffusion resistance of gaseous species in the pyrolysis process. 
For example, throughout the pyrolysis experiment of some typical 
Chinese coal, Chen et al. [48] and Yang et al. [49] observed that 
pressure had a great influence on coal pyrolysis and the solid char 
yield increased with increasing pressure. Pressurized condition 
would affect the chemical structure of char and further strengthen 
the volatile matter precipitation in the secondary reaction. Thus, 
the pressure effect should be accounted for in a realistic model. A 
mathematical model proposed by Anthony and Howard [50] cap¬ 
tures the pressure effect on the volatile products by considering 
the competition between diffusional escape and secondary reac¬ 
tion of reactive volatile species during a pyrolysis process. The 
model expression is given by: 

w*(P t ) = w* r + w; x (l + p t x k,/kc) (2) 

where w*, w* r and w** denote respectively the mass portion of vol¬ 
atile products at the studied pressure, very low pressure 
(<0.001 atm) and very high pressure (>100 atm), P t is the operating 
pressure (atm) and /ci//c c = 0.56. Wen and Chaung [51 ] modified this 
model and used a linear interpolation of the data of Anthony and 
Howard [50] to account for the pressure effect, as shown in the fol¬ 
lowing equation: 

w*(P t ) = w r * X (1 -ocxIn(P t )) (3) 

where w* is the mass portion of volatile products at ambient pres¬ 
sure and a is assumed to be approximately 0.066 for lignite and 
bituminous coal. 

In this work, a mathematical model that accounts for the pres¬ 
sure effect is established to describe the overall pyrolysis process 
in the Lurgi gasifier, as described below. 


reported that the dominant sulfur-containing product is 
H 2 S, and the remaining sulfur-containing product is COS 
[55]. H 2 S is the product of desulphurization reactions 
between H2, generated by the decomposition of the coal 
pyrolysis reaction, and pyrite sulfur [54,55]. Khan [56,57] 
investigated the pyrolysis of a number of coals at a relatively 
low temperature (770 K) in a fixed-bed reactor. The COS/H 2 S 
ratio was found to depend on the coal rank and sulfur con¬ 
tent, varying from 1:10 to 1:4. Considering that the lignite 
coal under study here is a low-rank coal with a high content 
of sulfur, the ratio of COS/H 2 S is assumed to be 1:4. Thus, it is 
assumed that 80% of the S atoms are consumed in the forma¬ 
tion of H 2 S. 

Wicoal FCam \ 

MWp^J ~ MWp^JJ 

16 

= "V x Na U + Y^ N Pyi x Na ‘j x ' = 1 2 * * * ( 6 ) 

J=1 j=2 

mu, = £> pyrj - x Na tj i <= [ 2 ,6] ( 7 ) 

Eqs. (8) and (10) model the distribution of elements N and S, 
respectively. The molar flow rate of NH 3 is identical to that of 
HCN, as shown in Eq. (9). 


OAmfcoal | 

w i;COd \ 

— YNpyrj x Nciij i — 3 

s=N 2 

(8) 

MWpyJ 

Npyr,NH 3 = 

Npyr,HCN 


( 9 ) 

O.Smfcoal | 

' w l)Coa/ \ 
,MW pyr J 

= N pyrj x Nciij i = 5 

s=H 2 S 

(10) 


(1) Firstly, the sets of all atoms At pyrti and vapor-phase products 
Vppyrj are defined as follows: 

e e At pyr j = {C, H, N, Cl, S, 0} i e [1,6] (4) 

5 c Vp pyrJ = {C s ,C0,C0 2 ,CH 4 ,C 2 H 6 ,H 2 ,N 2 ,H 2 S,C0S;NH 3 , 
HCN,C 9 H 10 , C 6 H 6 O,C 10 H 16 , C 9 H 12 , H 2 0}, j e [1,16] (5) 

where i and j denote each atom and gaseous species in the pyrolysis 

reaction, respectively. 

(2) The constraints of atom balance on C, N, Cl, S and O are 
shown in Eqs. (6) and (7). For the conversion of N atoms in 

the coal pyrolysis reaction, it has been reported that the 

major nitrogenous products are N 2 , NH 3 and HCN. The stud¬ 

ies of Baliban et al. [46] Liu et al. [52] and Di Nola et al. [53] 

indicate that 40% of N atoms in coal end up in N 2 and that 
HCN/NH 3 ^1. For the conversion of S atoms, it has been 


where m/, w, MW, Na and N represent the mass flow rate, mass frac¬ 
tion, molar weight, number of atoms and molar flow rate, respec¬ 
tively; MF am and w* are the mass fraction of fixed carbon (ash/ 
moisture free basis) and mass fraction based on experimental data, 
respectively. 

(3) The relationships between reference yields (see Table 5) and 
expected yields are treated as important constraints. Eqs. 
(11) and (12) illustrate the composition ratio constraints of 
coal gas and liquids while Eqs. (13) and (14) illustrate the 
constraints of distribution ratio relationships among coal 
gas, liquids and water. 


” mfcpyrl _ rfcpyrj ^ 
m fcpyr j+1 m fcpyrj+l 

m fcpyrj m fcpyr j ^ D 

- f -75- ^ Kfc 

, m JCpyrj +1 Oy C pyrj +1 


j e [ 2 , 5 ]; fee [ 1 , 4 ] 


( 11 ) 
















C. He et al./Applied Energy 2 2 2 (2013) 742-757 


749 


m /cpyiy _ m 2cpyrj ^ 

Wcpyrj+1 m Cpyrj +1 

m fcpyrj _ m fcpyrj < ^ 
mfcpyrj+1 mf C pyrj + l 


je [12,14]; z e [1,3] 


( 12 ) 


27i/; 


J2j=2 m fcpy 

y^j-19 Wlfcpyrj Cpyr,coal liquids 


Cpyr, coal gas 


» -G, 


X^j =2 m /cpyrj 
Ej!n m fcpyrj 


mf ; 
m/; 


Cpyr,coal gas 




Cpyr,coal liquids 


X)j=l 2 mfcpyrj Cpyr,coal liquids ^ 

Xyj=l 6 m fcpyr,j m fcpyr ;water 

J2j=n m fcpyrj _ m f Cpyr,coal liquids < 
J2j=!6 m fcpyr.j m fcpyr ,water 


(13) 


(14) 


r Rfe > 0, L z 3 s 0, 

l Ncpyrj > 0 


Gi2 ^ 0 


Vj c [1,16] 


(15) 


In Eqs. (11)—(14), we define three slack variables, called ft, L and G, 
which represent the constraints of composition ratio relationships 
between the expected pyrolysis products and reference pyrolysis 
products based on the experimental data listed in Tables 5 and 6 . 

(4) The objective function is given by 

4 3 2 

+ cpL z + yG r ) (16) 

k =1 z= 1 r=l 


In Eq. (16), the three weighting factors, 2 c, <p and 7 (k ^ 1, ^ 1, 

7 ^ 1 ), are used to denote the relative importance of the constraint 
coefficients ft, L and G, respectively. In the proposed model, 
k = (p = y = \. The pyrolysis model is optimized in the GAM 22.3 
modeling system. The resulting optimization problem has 42 vari¬ 
ables and 43 constraints with 209 Jacobian elements, 130 of which 
are non-linear. Solving this model to locate its global optimality by 
means of the global solver BARON gives the N pyrj values listed in Ta¬ 
ble 7. The final overall pyrolysis reaction for the given coal is ex¬ 
pressed in the following equation: 

C 6 . 197 H 4 . 524 N 0 . 062 S 0 . 069 O 1.123 —> 5.557C( S ) + 0.106CO 

+ O.O 8 OCO 2 + O.I 2 OCH 4 
+ 0.084C 2 H 6 + 0.419H 2 
+ 0.012N 2 + 0.019NH 3 
+ 0.014COS + 0.056H 2 S 
+ 0.019HCN + 0.008C 14 H 10 
+ 0.006C 6 H 6 0 + 0.006C 7 H 8 
+ O.O 2 OC 7 H 10 + 0.834H 2 O (17) 

The entire fixed carbon will be converted into char, and some 
distillable liquids will be decomposed in high pressure circum¬ 
stance, so that a total of 0.730 moles of carbon atoms in volatile 
matter is transferred to the char at the end of the pyrolysis 
reaction. 


3.2. Models of char gasification and combustion processes 

In the process of char gasification and combustion, three types 
of homogeneous reactions (gas-gas reactions) and four types of 
heterogeneous reactions (gas-solid reactions) are considered, see 
Table 2. The macro-reaction kinetics and reaction rate constants 
for these reactions are given in Table 8 . Among the four heteroge¬ 
neous reactions, reactions R2-R4 are volumetric reactions in which 


Table 7 

Fractional distribution of coal pyrolysis. 


Product 

Molar flow 

Product 

Molar flow 

FC a 

4.827 

C a 

0.730 

c 2 h 6 

0.084 

n 2 

0.012 

h 2 o 

0.834 

h 2 s 

0.056 

CO 

0.106 

cos 

0.014 

C0 2 

0.080 

C 14 H 10 X 10 3 

8.504 

ch 4 

0.120 

C 6 H 6 0 X 10 3 

6.118 

nh 3 

0.019 

C 7 H 8 X 10 3 

5.834 

HCN 

0.019 

C 7 H 10 X 10 3 

19.890 

h 2 

0.419 




a Fixed carbon (FC) and carbon (C) generated from pyrolysis are included in C (s) 
(char). 


gaseous matters can quickly diffuse into particles and the reactions 
take place throughout the interior of the particles. Reaction R5 by 
contrast is a surface reaction in which the overall reaction rate is 
influenced by the joint effect of the chemical reaction, diffusion 
in ash particles and gas film mass transfer. The unreacted shrinking 
core model is used to describe the reaction and mass transfer rates 
in reaction R5, as shown in the following equations: 


_ r 1 1 1 

S " Cliff + l< s Y 2 + Clash 

I} -1 )] 

(18) 

y = dcore =(1 Xch J 
Uparticle 


(19) 

kdash — kdiff ' &b 


( 20 ) 


where Y denotes the ratio of radius of unreacted core (d core ) in cm to 
average radius of the char particle (d par ticie) in cm; X char is the carbon 
conversion ratio in char (%); k dif f -, k s and k daS h respectively denote the 
gas film mass transfer coefficient, chemical reaction rate constant 
and ash diffusion coefficient (mol (cm 3 atm s) _1 ); P t and P- respec¬ 
tively represent the partial pressure of the reverse reaction and 
equilibrium pressure of gaseous species i (atm); s b denotes the 
porosity of ash (dimensionless). 

For countercurrent contact mode, a series of RCSTRs can be used 
to approximately simulate the mass transfer and chemical reaction 
of the gas-solid reactants [22]. The RCSTR has the characteristics 
that temperature and concentration are uniform and reactions pro¬ 
ceed instantaneously. Theoretically, the accuracy of the approxi¬ 
mation increases with increasing number of RCSTRs. However, 
iterative calculations based on a high number of CSTRs are usually 
time consuming and may exhibit slow solution convergence be¬ 
cause the form of the mathematical model of countercurrent 
fixed-bed gasifier is a two point boundary value problem [23]. As 
such, this work presents a novel solution strategy that considers 
the reactions of the gasification zone and combustion zone sepa¬ 
rately. This solution strategy can reduce the computation time 
and scale and the complexity of the algorithm for each RCSTR. Dif¬ 
ferent external FORTRAN subroutine codes have been developed to 
link m number of RCSTRs and n number of RCSTRs in Aspen Plus for 
simulating the gasification zone and combustion zone, respec¬ 
tively. Specifically, reactions R2-R4, R 6 and R7 are considered in 
the char gasification zone (from point A to point B in Fig. 3) while 
reactions R5, R7 and R 8 are considered in the char combustion 
zone (from point B to point C in Fig. 3). The interaction boundary 
of the gasification and combustion zones is located at point B 
where all the feedstock oxygen has been consumed, as shown in 
Fig. 3. Calculating the macro-reaction rate of the chemical reac¬ 
tions listed in Table 8 is the main task of the external FORTRAN 
subroutine codes. The programming method is available in the As¬ 
pen Plus User Models Guide [58]. Note that although these reaction 
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Table 8 

Reaction rate expressions and corresponding rate constants for the char gasification and combustion processes. 


Reaction no. 

Reaction rate 3 

Rate constant 3 

Units 

Ref. 

R2 b 

3000exp(=5S51)Cc(P H2 o-PH 2 o) 

n* _ PcoPh 2 

h 2 ° exp (17.29-1^12) 

mol (cm 3 s) -1 

[7,24] 

R3 b 

exp (-7.087 -5SP)C C (P H2 -P*h 2 ) 

p , _ r Pc« 4 I 0 - 5 

r H 2 0 - [exp(-13.43-30JJS)J 

mol (cm 3 s) -1 

[7,24] 

R4 b 

3000 exp (=sp)C c (Pco 2 -Pm 2 ) 

p* — Pco 

1 h 2° “ exp (20.92-2228Q) 

mol (cm 3 s) -1 

[7,24] 

R5 C 

[l.244(7/1800)* 75 + k dlr ef> (v 1 )] (P<J ’ > 

“ 

mol (cm 3 s) -1 

[7,24] 

R6 

312exp(^|m) 

- 

s -1 

[51] 

R7 

2.96 x 10- 4 ■ T 2 exp (- (c c0 ■ Ch 2 o - 

k W gs = exp (-4.3 + ^8) 

mol (m 3 s) -1 

[34,59] 

R8 

8.83 x 10 8 exp (=12003) (C Hz ) 2 

- 

mol (m 3 s) -1 

[34] 


a Q denotes the molar concentration (mol m -3 ) of reactant or product i; T denotes reaction temperature (K). 
b R2-R4 are slow chemical reactions and the overall reaction rate depends on chemical kinetics; kdiff ^ k s , kdash > k s (see Eq. (18)). 
c R5 is a fast chemical reaction and the overall reaction rate depends on mass transfer; k s ^ k d ijf, k s ^ k daS h (see Eq. (18)). 


rates and constants are applicable to the Lurgi fixed-bed coal gas¬ 
ifier under study, it is likely that testing of different kinetic data 
sets is required in order to select the most appropriate kinetic data 
for a given type of coal gasifier. 

Krishnudu et al. [60] divided a pressurized moving bed gasifica¬ 
tion pilot plant into 14 sequential layers along the height of the 
gasifier after analyzing the proximate and ultimate of the sequen¬ 
tial layer samples to represent three major zones: fast devolatiliza¬ 
tion, slow devolatilization and gasification-combustion. In 
practical situations, the height of the combustion-gasification 
zones accounts for about 80% of the total height of the gasifier 
and the height of the combustion zone at the bottom of the gasifier 
is usually shorter than that of the gasification zone, especially for 
young coal with higher reactivity (e.g., lignite and bituminous coal) 
which reacts rapidly with gasifying agents [45,61]. 

In Fig. 4 an iterative calculation procedure for determining the 
numbers of RCSTRs needed for the gasification and combustion 
zones is presented. The initial calculation condition assumes that 
each single RCSTR corresponds to a sequential layer with 0.45 m 
thickness along the height of the gasifier and is equivalent to a vol¬ 
ume of 5 m 3 . The first phase of the iterative calculation begins by 
increasing the number of RCSTRs in the combustion zone (n) and 
the gasification zone (m) simultaneously until the oxygen supply 
in the combustion zone is totally exhausted. Then the second phase 
of the iterative calculation takes over which continues to increase 
the number of RCSTRs in the gasification zone (m) until the amount 
of remaining char achieves a steady value. 

Fig. 5 shows the amounts of oxygen and char as functions of m 
and n for the conditions described below. Based on the Aspen Plus 
model formulation given in Table 4 and gasification conditions for 
Run 1 listed in Table 11, the iterative calculation results show that 
the amount of oxygen is practically zero when m = n = 7 (Fig. 5a, at 
point B). The optimal number of RCSTRs in the combustion zone is 
thus given by n = 7. At this point the overall char conversion is 
72.8% (Fig. 5b, at point B). With n fixed at 7, the char conversion in- 



Fig. 4. Iterative calculation procedure for determining RCSTRs numbers (m and n). 

creases with further increases in the number of RCSTRs (m) in the 
gasification zone, reaching a stable value of 98% when m = 13 
(Fig. 5b, at point C). When m = 13 and n = 7, it can be seen that 
new oxygen amount is also exactly exhausted when m = 7, 
although the entire oxygen amount curve (Fig. 5a) lies slightly be¬ 
low the curve for the case of m = n = 7. In summary, the optimum 
number of RCSTRs in the gasification zone and the combustion 
zone is respectively 13 and 7. 


r 


m reactors 


n reactors 

If 1 



A 


Gasification zone 


cJo . cJo 




0 2 /Air+Steam 
► Ash 


B 


Combustion zone 


C 


Fig. 3. Flowsheet diagram of RCSTRs in char gasification and combustion processes. 
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a. Number of RCSTRs (n) 



1 2 3 4 5 6 7 8 9 10 11 12 13 14 

b. Number of RCSTRs (m) 


Fig. 5. (a) Oxygen amount and (b) char amount in the gasification-combustion zones as functions of RCSTRs numbers (m and n). 


Table 9 

Temperature, pressure, exergy flow rate of the main streams of the coal gasifier and gasification system. 


No. 

Stream 

Input/output 

Temperature, K 

Pressure, atm 

Exergy flow rate 

1 

Coal 

Input 

298.15 

1 

15.01 MJkg -1 h -1 

2 

Syngas 

Output 

500 

40 

- 

3 

Steam 

Output 

450 

6 

0.76 MJ kg” 1 h” 1 

4 

Electricity 3 

Input 

- 

- 

0.003 kW kg coal- 1 h 1 

5 

Ash + uc 

Output 


- 

0 

6 

Water 

Input 

298.15 

1 

0 

7 

Steam 

Output/input 

550 

40 

1.15 MJkg” 1 h” 1 

8 

Water 

Input 

298.15 

1 

0 

9 

Heat 

Input 

- 

- 

2.50 MJ kg water -1 h 1 

10 

Air 

Input 

298.15 

1 

0 

11 

Electricity 

Input 

- 

- 

0.42 kW N m 3 oxygen -1 h -1 

12 

Steam 

Input 

800 

90 

1.50 MJkg -1 h -1 

13 

Heat 

Input 

- 

- 

0.012 MJ kg oxygen -1 h -1 


a Electricity consumption of gas compressor (coal locker), coal distributor, ash distributor, etc. 


3.3. Heat transfer model 

The heat transfer (Qt r ) from the gasifier to the exterior wall is 
usually about 0-30% of the heat value of the coal fed to the gasifier 
and is mainly due to heat transfer from the combustion zone wall 
(Q.t r om ), heat transfer from the gasification zone wall (Qf r asi ) and heat 
transfer from the drying-pyrolysis zone wall (Q^), as indicated in 
the following equations: 


Q-tr — 

or + Qf' + Q? 

(21) 

n c°m 

"tr 

m 

*-Uj2MTi-T wall ) 

i= 1 

(22) 

O gasi ; 

* - Twall) 

i= 1 

(23) 

o dp * 

^tr 

* Adp(Tdp - Twall) 

(24) 


where U denotes average heat transfer coefficient (J(m 2 hl<) *) 
which is assigned a value of 170 J (m 2 h K) -1 ; A is axial area of each 


Table 10 

Operating conditions used for simulating the pressurized Lurgi gasifier. 


Parameter 

Run 1 

Run 2 

Run 3 

Run 4 

Ambient conditions 

T = 298 K 

T = 298 K 

T= 298 K 

T = 298 K 


P = 1 atm 

P = 1 atm 

P = 1 atm 

P = 1 atm 

Bed size 

h = 12.5 m 

h = 12.5 m 

h = 12.5 m 

h = 12.5 m 


d = 3.8 m 

d = 3.8 m 

d = 3.8 m 

d = 3.8 m 

Feed coal 

40.28 t h -1 

40.28 t h -1 

36.65 th -1 

36.65 th -1 


T = 298 I< 

T = 298 K 

T = 298 K 

T = 298 K 

Feed steam 

S/C a = 0.795 

S/C = 0.752 

S/C = 0.797 

S/C = 0.752 


T= 708 I< 

T= 708 K 

T= 708 K 

T= 708 K 


P = 40 atm 

P = 40 atm 

P = 35 atm 

P = 35 atm 

Feed 0 2 (95.8 vol.%) 

0/ 

0/ 

0/ 

0/ 


C b = 0.186 

C = 0.179 

C = 0.188 

C = 0.179 


T= 383 K 

T = 383 K 

T= 383 K 

T = 383 K 


P = 40 atm 

P = 40 atm 

P = 35 atm 

P = 35 atm 

Operating pressure 

40 atm 

40 atm 

35 atm 

35 atm 

Wall temperature 

530 K 

530 K 

530 K 

530 K 


a S/C: steam to coal mass ratio. 
b O/C: oxygen to coal mass ratio. 
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Table 11 

Comparison of model predictions (Cal.) and industrial data (Ind.). 


Vapor comp. 

(yO 

Run 1 


Run 2 


Run 3 


Run 4 


Basic index (jb,) 

Run 1 


Run 2 


Run 3 


Run 4 


Cal. 

Ind. 

Cal. 

Ind. 

Cal. 

Ind. 

Cal. 

Ind. 

Cal. 

Ind. 

Cal. 

Ind. 

Cal. 

Ind. 

Cal. 

Ind. 

CO 

15.05 

14.47 

14.88 

15.42 

15.21 

14.94 

15.03 

15.44 

LHVg, MJ N ITT 3 

11.83 

11.98 

11.79 

11.55 

11.57 

11.69 

11.63 

11.68 

C0 2 

32.6 

32.15 

32.55 

33.63 

32.94 

32.51 

32.76 

32.52 

RE C ,% 

1.25 


2.08 


1.03 


0.43 


H 2 

38.45 

39.26 

38.93 

37.63 

38.53 

39.14 

38.94 

38.71 

4cg,% 

79.33 

79.24 

81.04 

78.44 

76.24 

78.01 

79.26 

76.74 

ch 4 

12.53 

12.84 

12.35 

11.98 

11.8 

12.01 

11.88 

11.95 

RE C ,% 

0.11 


3.31 


2.27 


3.28 


nh 3 

0.18 

0.17 

0.17 

0.19 

0.22 

0.19 

0.21 

0.2 

V g , 1<N m 3 kg” 1 

1.28 

1.25 

1.28 

1.27 

1.24 

1.23 

1.25 

1.24 

HCN 

0.18 

0.17 

0.17 

0.19 

0.22 

0.19 

0.21 

0.2 

RE C ,% 

2.40 


0.79 


0.81 


0.80 


c 2 h 6 

0.11 

0.1 

0.1 

0.1 

0.13 

0.11 

0.12 

0.11 

X* % 

99.70 

98.60 

99.50 

98.30 

98.00 

98.20 

99.30 

97.80 

n 2 

0.15 

0.14 

0.14 

0.14 

0.16 

0.16 

0.14 

0.15 

RE c ,% 

1.12 


1.22 


0.20 


1.53 


h 2 s 

0.57 

0.55 

0.54 

0.56 

0.62 

0.57 

0.54 

0.56 

m Ash , kg kg -1 

0.23 

0.25 

0.24 

0.26 

0.25 

0.27 

0.26 

0.29 

cos 

0.18 

0.15 

0.16 

0.16 

0.17 

0.18 

0.17 

0.16 

RE C ,% 

8.00 


7.69 


7.41 


10.35 


Tar 3 , g N m -3 

9.13 

8.34 

9.13 

8.34 

9.53 

8.62 

9.53 

8.62 

Tsyngasi K 

578 

508 

568 

508 

553 

498 

553 

492 

Phenol 3 , g N m -3 

4.12 

3.93 

4.12 

3.93 

4.32 

3.97 

4.32 

3.97 

RE C , % 

13.78 


11.81 


11.04 


12.40 


Naphtha 3 , g N m -3 

3.84 

3.54 

3.84 

3.54 

3.98 

3.65 

3.98 

3.65 

Tpeaky K 

1362 

1230 

1350 

1204 

1327 

1189 

1307 

1170 

Oil 3 , g N m -3 

13.38 

12.68 

13.38 

11.68 

13.58 

12.79 

13.58 

12.79 

RE C , % 

10.73 


12.13 


11.61 


11.71 


ME h , % 

7.91 


6.69 


9.47 


6.24 


Qtr, MJ kg coal -1 

0.22 

- 

0.22 

- 

0.20 

- 

0.18 

- 










Qioss, MJ kg coal -1 

0.19 

- 

0.19 

- 

0.178 

- 

0.16 

- 










Qjacket, kj kg coal -1 

26.8 

- 

25.3 

- 

24.2 

- 

22.6 

- 


a Industrial data of high-weight hydrocarbons (tar, phenol, naphtha and oil) vary with the mass flow of by-products in downstream units (gas water separation unit, phenol 
recovery unit, Rectisol unit). 

b me = v / s!l,[(yL <f -yL II )/J'L] 2 /N [66]. 
c RE»lbi„ d -b^/4f- 



Fig. 6. Exergetic analysis model of the coal gasifier process and gasification system. 


zone (m 2 ); T it T dp and T wa u respectively denote temperature of the 
2 th RCSTR in the gasification zone and combustion zone, average 
temperature of the drying-pyrolysis zone and average temperature 
of exterior wall (K). T dp and T wa u are assumed to be 700 K and 530 K, 
respectively. 

The heat is transferred to cooling water in the water jacket 
(Qjacket ) and to the surroundings as heat loss (QiossY The heat trans¬ 
ferred to the water jacket may be used to generate part of the 
steam needed for gasification. According to industrial data, the 
quantity of medium pressure steam generated from the water 
jacket is about 10% of the total gasifying steam (Qg S ). A heat balance 
gives the heat loss, as shown in the following equation: 

Qloss — Qtr ~ Qjacket = Qtr ~ O lQgs (25) 

Additional heat may be recovered from the crude syngas which 
is sent to the scrubbing cooler and cooled down to 473 K. The 
cooled syngas is then sent to the heat recovery steam generator 
(HRSG) to produce low pressure steam (6 atm), as shown in Fig. 2. 


4. Thermodynamic performance indicators 

The conversion efficiency of coal to syngas can be defined in 
numerous ways, including (1) feedstock conversion efficiency; 
(2) syngas yield efficiency; and (3) energy-based or exergy- 
based efficiency. Among the above indexes, the cold gas effi¬ 
ciency {vjcg) is considered a vital index for accessing gasifier 
performance which can be calculated from the following 
formula: 

_ syngas ~ LHV syng as 
Wlffeed ' LHVf eed 

where mf S y ngas and mff eed are the mass flow rate of syngas and feed 
fuels (kg h -1 ), respectively. 

However, the cold gas efficiency ( rj C c ) does not consider the ef¬ 
fect of process streams (steam and oxygen) and electricity con¬ 
sumption, and the effect of the waste heat recyclability of crude 
syngas is also neglected. According to the studies of Prins et al. 
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Fig. 7. Predicted vapor composition results with error bars from the industrial data. 
For a clearer illustration, the compositions of NH 3 , HCN, C 2 H 6 , N 2 , H 2 S and COS have 
been multiplied by a factor of 10. Other results (e.g., LHV g , vj C g and m Ash ) are not 
shown (see Table 11). 


[29,62,63] the exergetic efficiency is recommended for assessing 
the material/energy conversion efficiency of solid feedstock (coal, 
biomass, solid waste, etc.) gasification. 

4.2. Exergetic efficiency in coal gasification 

For conventional coal gasifiers, most researchers use the ratio of 
the exergy of outgoing products to that of incoming fuels to calcu¬ 
late the exergetic efficiency /71 [31], as shown in the following 
equation: 


m = 


E syngas(E XS y n g as + syngas) 


m ffuel(E x j ue i) + mfsteam(E C xsteam + ^x,steam) + m foxygen (E x ,oxygen + ^x,o 


'x.oxyge 


(27) 


where Ef xfyngas ,E% eam ,E c f oxygen , Keygen and E f S uei represent the chemi- 
cal exergy of crude syngas, steam, oxygen and the feedstock coal, 
respectively; Ef syngas ,E^ leam and E p x h oxygen represent the physical exer- 
gy of crude syngas, steam and oxygen, respectively. 

Although the unreacted char in the gasifier is not a target prod¬ 
uct, it has both chemical and physical exergy in each reaction zone 
and should be included in the calculation of the exergetic effi¬ 
ciency. For example, the unreacted char leaving the gasification 
zone is a useful outgoing flow which may be used as feedstock in 
the subsequent combustion zone and converted to syngas via par¬ 
tial oxidation reaction. Additionally, the exergetic efficiency should 
also consider the sensible heat transferred from gasifier wall to 
heat cooling water in the water jacket and generate considerable 
amounts of medium pressure steam (550 K, 40 atm) which may 
be used as gasifying medium. As such, Eq. (27) is modified to in¬ 
clude the exergy of unreacted char mf uc (E x h uc + E ph uc ) and jacket 
steam mf jacke[ (E c x h saam + E p x h steam ) in the output exergy, as shown in 
the following equation: 


% = 


E x ,gas_out E X g as j n 
Ex.fuelJn — E x f ue i_ uc 


E m f syngas (E ; 


’ch 


syngas \^x,syngas ' 


" ^x,syngas) + m fjacket (Exsteam ' 


-E ph ) 

x,steam J 




m fsteam (E XjS team 


' Ex,steam) + m foxy gen (E Xfi 


+ E ph ) 

,oxygen 1 x,oxygen) 


mfU^ue,)-rnfuc(C c +E^ c ) 


(29) 


However, the coal gasifier is just one part of a gasification sys¬ 
tem, albeit a major one. As shown in Fig. 6 , the complete gasifica¬ 
tion system also consists of auxiliary facilities such as the air 
separation unit (ASU), heat recovery steam generation unit (HRSG), 
water jacket, steam generator, etc. Irreversibility occurs in these 
units since they consume or recover exergy. The overall exergetic 
efficiency is therefore given by: 


4 overall 


E m f syngas (E x ,syngas ""I - 


Ex,syngas ) ( E f,steam ' 




m ffuel(E x f ue i) + E W electricity 


(30) 


where I x , steam and \ Xt0 xy g en (kj kg coal -1 ) are the irreversibility (inter¬ 
nal exergy loss) in the steam generators (Heating 1 and Mixing 1 
and 2 in Fig. 6 ) and oxygen generators (ASU and Heating 2 in 
Fig. 6 ), respectively. W e i ectricit y (kw kg coal -1 ) is power exergy for 
normal operation of the gasifier (such as coal distributor and ash 
distributor). 

To determine the chemical exergy of raw coal, char and fuel gas, 
a literature correlation [65] is used: 


E c x h = px LHV 


(31) 


'1.0438 + 0.0158x|ffl + 0.0813x^l, < 0.5 

/’ = < 1.0414+0.0177x!g-0.3328x5gx(l+0.0537x!^ 

1 -0.4021 x|g 

(32) 

where LHV is the lower heating value of fuel which is given by 
LHV = 0.943 x HHV where HHV is the higher heating value. HHV 
(volume basis) can be calculated by the HHV- 0 function using the 
REFSYS correlation property from the Aspen Plus Physical Property 
Database; n(C), n(H) and n(O) are the mole fractions of carbon, 
hydrogen and oxygen in the fuel, respectively. Additionally, the 
chemical exergy of each high-weight hydrocarbon is assumed to 
be equal to that of benzene (C 6 H 6 ). 

The physical exergy of pure substance is given by: 

= (H- H 0 ) - T 0 (S - So) = (H- T 0 S) - (H 0 - T 0 S 0 ) (33) 

where H and S are enthalpy and entropy at given temperature and 
pressure, and H 0 and S 0 are the values of these parameters at ambi¬ 
ent temperature and pressure (298.15 K, 1 atm). The thermody¬ 
namic availability for a mixture, H-T 0 S and H 0 -T 0 S 0 , can be 
calculated from the Aspen Plus Physical Property Database using 
the AVAILMX function (T 0 = 298.15 K). 

The temperature, pressure and exergy flow rate of the main 
streams of the coal gasifier and gasification system are listed in 
Table 9. 


0 5 < ”(2) < 2 

U J ^ n(C) ^ z 


^2 


m f syngas (E x h syngas 


™ffuel(E C x [ 


'x fuels 


+ ^ amgas ) + mfuc{^ 

mfsteam (E xf team 


rph 

^x, steam. 


-p + mWEg 

+ mf oxygen {E Xl0 xygen 


-ch 

,steam " 


rph 

ctP 


-E ph ) 

x,oxygen J 


( 28 ) 


According to the characteristics of fossil fuel conversion, Grass- 
man proposed another exergy-based index which is known as the 
intrinsic or rational efficiency [64]. As expressed in Eq. (29), it is 
defined as the exergy increase of gas divided by the exergy de¬ 
crease of solid fuel. 


5. Results and discussion 

5.1. Model validation 

The model is validated using data obtained from the industrial- 
scale Lurgi gasifier. The NLP model’s results and reaction kinetics 
have been inserted into the simulation, and input parameters of 
the Aspen Plus model are listed in Table 3. Besides, the operation 
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conditions for the pressurized Lurgi gasifier are listed in Table 10. 
Run 1 is the steady state operating data of this SNG plant and Runs 
2-4 are basic plant design data. 

Fourteen syngas component concentrations, ten performance 
indexes and error analysis based on the reference data are listed 
in Table 11. The model predictions are in good agreement with 
industrial measurements for gas component concentrations (y f ), 
carbon conversion efficiency (X c ), low heating value ( LHV ), dry 
gas yield (vy, and cold gas efficiency (2 /cg)- The mean error (ME) 


is below 10% and the largest relative error (RE) is less than 15%. 
Large relative errors (RE) are observed for both the predicted peak 
bed temperatures and outgoing syngas temperatures which are 
about 20-100 K higher than the industrial data. The discrepancies 
may be attributed to imperfect gas-solid mixing in the Lurgi gas¬ 
ifier which would violate the model’s assumption of perfect gas- 
solid mixing. Based on the results in Table 11, the predicted vapor 
composition results with error bars from the industrial data are 
presented in Fig. 7. Using Run 1 as example, Fig. 8 shows the con- 



Fig. 8. Profiles of gas component/char concentrations and temperature as a function of gasifier bed height for Run 1. 


(a) Combustion zone 


(c) Lurgi gasifier 



Oxygen/Coal ratio (kg/kg) 



Oxygen/Coal ratio (kg/kg) 


(b) Gasification zone 
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<D 
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30.0 


20.0 

0.15 0.165 0.18 0.195 0.21 0.225 

Oxygen/Coal ratio (kg/kg) 



(d) Gasification system 



Fig. 9. Effects of oxygen/coal ratio on the exergetic efficiency and exergy loss of (a) combustion zone, (b) gasification zone, (c) Lurgi gasifier, and (d) entire gasification system 
(feed coal: 40.28 t h _1 ). 
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(a) Combustion zone 



Steam/Coal ratio (kg/kg) 




w 


(b) Gasification zone 



Steam/Coal ratio (kg/kg) 


(c) Lurgi gasifier 



Steam/Coal ratio (kg/kg) 


(d) Gasification system 



Steam/Coal ratio (kg/kg) 


Fig. 10. Effects of steam/coal ratio on the exergetic efficiency and exergy loss of (a) combustion zone, (b) gasification zone, (c) Lurgi gasifier, and (d) entire gasification system 
(feed coal: 40.28 t h -1 )- 


Table 12 

Thermodynamic performance of the gasification system. 



Case 1 

Case 2 

Case 3 

Pressure, atm 

40 

35 

30 

Overall exergetic efficiency, % 

72.4 

70.1 

65.2 

Total exergy loss, MJ kg coal -1 

5.67 

5.88 

6.16 


centrations of some gas components and char and temperature as a 
function of bed height. The much lower concentrations of high- 
weight hydrocarbons are not shown in Fig. 8 . 

5.2. Thermodynamic analysis 

After the model is validated with industrial data, it is used to 
investigate the effects of oxygen/coal ratio and steam/coal ratio 
on the gasification performance. In the simulations, only the oxy¬ 
gen/coal or steam/coal ratio is varied; all the remaining input 
parameters retain the values given for Run 1 (see Table 10). 

The oxygen/coal ratio is a key parameter that controls the oper¬ 
ating temperature and carbon conversion rate in the Lurgi gasifier. 
Fig. 8 illustrates different exergetic efficiency profiles such as r\ A , 
r\ 2 , 273 , and rj overall calculated according to Eqs. (27)-(30) as well 
as exergy loss profiles at different oxygen/coal ratios for the com¬ 
bustion zone, the gasification zone, the Lurgi gasifier and the over¬ 
all gasification system. Comparison of Fig. 9a and b indicates that 
increasing the oxygen/coal ratio has a more significant impact on 
the exergetic efficiency of the combustion zone than that of the 
gasification zone. 


When the oxygen/coal ratio increases from 0.150 to 
0.225 kg kg -1 , the exergy loss profile of the combustion zone goes 
through a minimum (Fig. 9a). Nevertheless, the exergy loss is sig¬ 
nificantly higher than that of the gasification zone which is in 
the 1.10-1.30 MJ kg -1 range (Fig. 9b). Fig. 9d shows that the exergy 
loss profile of the overall gasification system also goes through a 
minimum, varying from 7.30 MJ kg -1 to 8.20 MJ kg -1 . Fig. 9 illus¬ 
trates that the highest exergy loss occurs in the combustion zone 
which accounts for over 60% and 40% of the exergy loss of the Lurgi 
gasifier and the overall gasification system, respectively. It is inter¬ 
esting to observe in Fig. 9b that the exergetic efficiency r\ A of the 
gasification zone exhibits relatively low values (48-55%), which 
are smaller than those of the combustion zone and the overall gas¬ 
ification system. This is due to the fact that the exergy available in 
the unreacted char is not considered by the indicator 271 . This unre- 



Fig. 11. Exergy loss of different units within the gasification system. 
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acted char flows to the combustion zone but does not undergo fur¬ 
ther reaction. The unreacted char is ultimately discharged from the 
gasifier. The exergetic efficiencies 271 and rj 2 of the combustion zone 
and the overall gasification system are therefore nearly the same 
(see Fig. 9a and c). 

In Fig. 9d, it can be seen that the overall exergetic efficiency of 
the entire gasification system ( f] overa ii ) increases and then decreases 
with increasing oxygen/coal ratio. The maximum rj ovemU is found at 
an oxygen/coal ratio of about 0.195. The initial increase in y] overa u 
with increasing oxygen/coal ratio is due to the enhancement of 
char conversion in the combustion zone by the presence of more 
oxygen. Additionally, the heat generated from the char oxidation 
will facilitate the endothermic reduction reactions occurring in 
the gasification zone, boosting further the value of 2 / overa //. 

Once the oxygen/coal ratio goes beyond a certain value, the CO 
combustion (CO + 0 2 -► C0 2 ) in reaction R5 is enhanced, causing 
the concentration of noncombustible C0 2 to increase and the con¬ 
centrations of other combustible gas species (e.g., hydrogen and 
methane) to decrease. Consequently, r\ ove raii decreases due to the 
fact that the decrease in the chemical exergy of the vapor-phase 
products is not fully compensated by the increase in their physical 
exergy. Additionally, increasing the oxygen feedstock increases 
electricity consumption in the air separation unit (ASU), lowering 
^1 overall- Because uncertain fluctuations of oxygen feed may occur 
in practical operations, the oxygen/coal ratio should be kept within 
the range of 0.180-0.190, and it must be not higher than 0.195. 
This is because excess oxygen feed may lead to an abnormal oper¬ 
ation temperature and even a serious slagging problem in the gas¬ 
ifier bed (T>1500K, see Table 3). In summary, the optimum 
oxygen/coal ratio is found to be in the region of 0.19. 

Using the optimum oxygen/coal ratio of 0.19 and a coal feed of 
40.28 th -1 , additional simulations are conducted by varying the 
steam/coal ratio from 0.45 to 1.25. The results are shown in 
Fig. 10. The effect of steam/coal ratio on the exergetic efficiency 
of the gasification zone (Fig. 10b) is more pronounced than that 
of the combustion zone (Fig. 10a). As can be seen in Fig. 10b, the 
three exergetic efficiencies 271 , ri 2 and 273 first increase and then de¬ 
crease with increasing steam/coal ratio. Initial increases in the 
steam/coal ratio enhance the char gasification (C + H 2 0 -> CO + H 2 ) 
and CO shift (CO + H 2 0 -> C0 2 + H 2 ) reactions, increasing the yields 
of gas components (methane, carbon monoxide and hydrogen) in 
the gasification zone which result in higher exergetic efficiency. 
However, when the char gasification reaction reaches equilibrium, 
further increases in the oxygen/coal ratio exert little influence on 
the syngas production in the gasification zone. The reduction in 
exergetic efficiency at high steam/coal ratios is due to the fact that 
high steam concentrations reduce the operating temperature, ad¬ 
versely affecting performance of the gasification zone. For the Lurgi 
gasifier, Fig. 10c reveals that all three exergetic efficiency indica¬ 
tors reach maximum values at a steam/coal ratio of 0.80. The over¬ 
all exergetic efficiency of the entire gasification system also 
reaches a maximum value at the same steam/coal ratio, as shown 
in Fig. lOd. It should be noted that in practical operations the 
steam/coal ratio is usually between 0.80 and 0.90 in order to en¬ 
sure high carbon conversion rates. Hence, under the given condi¬ 
tions of Run 1, the optimum oxygen/coal ratio and steam/coal 
ratio are found to be 0.19 and 0.85, respectively. At these optimum 
ratios, the overall exergetic efficiency of the entire gasification sys¬ 
tem is slightly above 70%. 

Apart from the oxygen/coal ratio and the steam/coal ratio, the 
operating pressure of the Lurgi gasifier is also a key parameter 
affecting the thermodynamic efficiency of the gasification system. 
Three cases are simulated to investigate the influence of operating 
pressure at the gasification conditions of Run 1 and optimum oxy¬ 
gen/coal ratio and steam/coal ratio of 0.19 and 0.85, respectively. 
As can be seen in Table 12, decreasing the operating pressure from 


40 to 30 atm reduces the overall exergetic efficiency of the gasifi¬ 
cation system by about 7%. In addition, the exergy loss increases 
with decreasing operating pressure. Fig. 11 illustrates the exergy 
loss of different components within the gasification system for 
the three cases. It is evident that the higher exergy loss in case 3 
is attributable to increases in the exergy loss of the combustion 
and gasification zones. 

6. Conclusions 

We have developed a steady state kinetic model for an indus¬ 
trial-scale pressurized Lurgi gasifier using the simulator Aspen 
Plus. The thermodynamic efficiencies of the gasifier and the entire 
gasification system have been investigated in considerable detail 
via simulation studies. Five sequential modules: drying zone, pyro¬ 
lysis zone, gasification zone, combustion zone and overall heat 
recovery unit are considered. A stoichiometry-based mathematical 
model of the pyrolysis process and an external FORTRAN subrou¬ 
tine that predicts the syngas composition of the Lurgi gasifier are 
developed and integrated into the process simulation model. Satis¬ 
factory agreement between predicted results and industrial data is 
obtained, demonstrating the power of the proposed gasifier model. 

The model is used to investigate the effects of two key operating 
parameters, namely oxygen/coal mass ratio and steam/coal mass 
ratio, on the exergetic efficiencies of the gasifier and the gasifica¬ 
tion system as a whole. The simulation results show that the oxy¬ 
gen/coal ratio exerts a strong impact on the exergetic efficiency 
profile of the combustion zone within the Lurgi gasifier. By con¬ 
trast, the exergetic efficiency profile of the gasification zone is 
somewhat insensitive to the oxygen/coal ratio. The opposite trend 
is however observed in the case of the steam/coal ratio. The exer¬ 
getic efficiency of the gasification zone is more sensitive to the 
steam/coal ratio compared to that of the combustion zone. The 
overall exergetic efficiency of the entire gasification system 
reaches a maximum value of about 72.4% at an oxygen/coal ratio 
of 0.19 and a steam/coal ratio of 0.85. 
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